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Parity violating elastic electron scattering on neutron-rich nuclei is de- 
scribed in the framework of relativistic mean-field theory. Self-consistent 



■ ground state density distributions of Ne, Na, Ni and Sn isotopes are cal- 

culated with the relativistic Hartree- Bogoliubov model, and the resulting 



0^ ' neutron radii are compared with available experimental data. For the elastic 

■ scattering of 850 MeV electrons on these nuclei, the parity- violating asymme- 

Q . try parameters are calculated using a relativistic optical model with inclusion 

^ . of Coulomb distortion effects. The asymmetry parameters for chains of iso- 



topes are compared, and their relation to the Fourier transforms of neutron 



^ . densities is studied. It is shown that parity- violating asymmetries are sensi- 

tive not only to the formation of the neutron skin, but also to the shell effects 
of the neutron density distribution. 
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I. INTRODUCTION 



Measurements of ground state distributions of nucleons provide fundamental nuclear 
structure informations. The ground state densities reflect the basic properties of effective 
nuclear forces, and an accurate description of these distributions presents the primary goal 
of nuclear structure models. 

Extremely accurate data on charge densities, and therefore on proton distributions in 
nuclei, are obtained from elastic scattering of electrons. Experimental data of compara- 
ble precision on neutron density distributions are, however, not available. It is much more 
difficult to measure the distribution of neutrons, though more recently accurate data on dif- 
ferences in radii of the neutron and proton density distributions have been obtained. Various 
experimental methods have been used, or suggested, for the determination of the neutron 
density in nuclei [|l| . Among them, one that is also very interesting from the theoretical point 
of view, is parity violating elastic electron scattering. In principle, the elastic scattering of 
longitudinally polarized electrons provides a direct and very accurate measurement of the 
neutron distribution. It has been shown that the parity-violating asymmetry parameter, 
defined as the difference between cross sections for the scattering of right- and left-handed 
longitudinally polarized electrons, produces direct information on the Fourier transform of 
the neutron density 0. 

In a recent article |^, Horowitz has used a relativistic optical model to calculate the 
parity-violating asymmetry parameters for the elastic scattering of 850 MeV electrons on a 
number of spherical, doubly closed-shell nuclei. Ground state densities were calculated using 
a three-parameter Fermi formula and the relativistic mean-field model. Coulomb distortion 
corrections to the parity-violating asymmetry were calculated exactly. It has been shown 
that a parity violation experiment to measure the neutron density in a heavy nucleus is 
feasible. 

Informations about the distribution of neutrons in nuclei should also constrain the isovec- 
tor channel of the nuclear matter energy functional. A correct parameterization of the 
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isovector channel of effective nuclear forces is essential for the description of unique phe- 
nomena in exotic nuclei with extreme isospin values. For neutron-rich nuclei, these include 
the occurrence of nuclei with diffuse neutron densities, the formation of the neutron skin 
and the neutron halo. At the neutron drip-line different proton and neutron quadrupole 
deformations are expected, and these will give rise to low-energy isovector modes. 

In the present work we describe parity violating elastic electron scattering on neutron- 
rich nuclei in the framework of relativistic mean-field theory. The ground state density 
distributions will be calculated with the relativistic Hartree-Bogoliubov (RHB) model. This 
model represents a relativistic extension of the Hartree-Fock-Bogoliubov (HFB) framework, 
and it has been successfully applied in studies of the neutron halo phenomenon in light nuclei 
01, properties of light nuclei near the neutron-drip line ground state properties of Ni 
and Sn isotopes [Q, the deformation and shape coexistence phenomena that result from the 
suppression of the spherical N=28 shell gap in neutron-rich nuclei H, the structure of proton- 



rich nuclei and the phenomenon of ground state proton emission p|-p!0|. In particular, it has 
been shown that neutron radii, calculated with the RHB model, are in excellent agreement 
with experimental data In studies of the structure of nuclei far from the /3-stability 

line, both on the proton- and neutron-rich sides, it is important to use models that include a 
unified description of mean- field and pairing correlations [0. For the ground state density 



distributions of neutron-rich nuclei in particular, the neutron Fermi level is found close to the 
particle continuum and the lowest particle-particle modes are embedded in the continuum. 
An accurate description of neutron densities can only be obtained with a correct treatment 
of the scattering of neutron pairs from bound states into the positive energy continuum. 

Starting from the RHB self-consistent ground state neutron densities in isotope chains 
that also include neutron-rich nuclei, we will calculate the parity-violating asymmetry pa- 
rameters for the elastic scattering of 850 MeV electrons. The main point in the present 
analysis will be to determine how sensitive the asymmetry parameters are to the varia- 
tions of the neutron density distribution along an isotope chain. An interesting question is 
whether parity violating electron scattering could be used, in principle at least, to measure 



the neutron skin, or even the formation of the neutron halo. Of course, studies of electron 
scattering from exotic nuclei require very complex experimental facilities, as for example the 
double storage ring MUSES, under construction at RIKEN. By injecting an electron beam, 
generated by a linear accelerator, in one of the storage rings, and storing a beam of unstable 
nuclei in the other, collision experiments between radioactive beams and electrons could be 
performed. 

In Section II we present an outline of the relativistic Hartree-Bogoliubov model and 
calculate the self-consistent ground state neutron densities of Ne, Na, Ni and Sn isotopes. 
For the elastic scattering of longitudinally polarized electrons on these nuclei, in section III 
we use a relativistic optical model to calculate the parity-violating asymmetry parameters. 
Coulomb distortion corrections are included in the calculation, and the resulting asymmetries 
are related to the Fourier transforms of the neutron densities. The results are summarized 
in Section IV. 

II. RELATIVISTIC HARTREE-BOGOLIUBOV DESCRIPTION OF GROUND 

STATE DENSITIES 

In the framework of the relativistic Hartree-Bogoliubov model |]8|Jl^, the ground state 
of a nucleus |$ > is represented by the product of independent single-quasiparticle states. 
These states are eigenvectors of the generalized single- nucleon Hamiltonian which contains 
two average potentials: the self-consistent mean-field F which encloses all the long range 
particle-hole (ph) correlations, and a pairing field A which sums up the particle-particle 
(pp) correlations. The single-quasiparticle equations result from the variation of the energy 
functional with respect to the hermitian density matrix p and the antisymmetric pairing 
tensor k. The relativistic Hartree-Bogoliubov equations read 

hD-m~~X A \ /Uk{r)\ /Uk{r)\ 

-A* -ho + m + x) \Vk{r)) \Vkir) ) ' 
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hi) is the single-nucleon Dirac Hamiltonian, m is the nucleoli mass, and A denotes the pairing 
field. The column vectors represent the quasi-particle spinors and Ek are the quasi-particle 
energies. The chemical potential A has to be determined by the particle number subsidiary 
condition in order that the expectation value of the particle number operator in the ground 
state equals the number of nucleons. 



The Hamiltonian he describes the dynamics of the relativistic mean-field model JTS 
nucleons are described as point particles; the theory is fully Lorentz invariant; the nucleons 
move independently in the mean fields which originate from the nucleon-nucleon interaction. 
Conditions of causality and Lorentz invariance impose that the interaction is mediated by 
the exchange of point-like effective mesons, which couple to the nucleons at local vertices: 
the isoscalar scalar cr-meson, the isoscalar vector tu-meson, and the isovector vector p-meson. 
The single-nucleon Hamiltonian in the Hartree approximation reads 

(1 — T-) 

hD = -ia-V + I3{m + g^a{Y)) + ^^cu(r) + gpT^p-i{r) + ^ A{y), (2) 

where cr, tu, and p denote the mean-field meson potentials, g^, g^ and gp are the correspond- 
ing meson-nucleon coupling constants, and the photon field A accounts for the electromag- 
netic interaction. The meson potentials are determined self-consistently by the solutions of 
the corresponding Klein-Gordon equations. The source terms for these equations are sums 
of bilinear products of baryon amplitudes, calculated in the no-sea approximation. 
The pairing field A in the RHB single-quasiparticle equations (|l]) is defined 

A„b(r, r) = ^ XI K6cd(r, r')Kcd(r, r'), (3) 

where a, 6, c, d denote quantum numbers that specify the Dirac indices of the spinors, 
Vabcdi'^^i r') are matrix elements of a general two-body pairing interaction, and the pairing 
tensor is 

«:,,(r,r')= E U:,{v)V,,{v') . (4) 

Ek>0 

The input parameters of the RHB model are the coupling constants and the masses for 
the effective mean-field Hamiltonian, and the effective interaction in the pairing channel. In 



most applications we have used the NL3 set of meson masses and meson-nucleon couphng 
constants |0 for the effective interaction in the particle-hole channel: m = 939 MeV, 
m„ = 508.194 MeV, = 782.501 MeV, nip = 763.0 MeV, = 10.217, g2 = -10.431 
fm"^, = —28.885, Quj = 12.868 and Qp = 4.474. Results of NL3 model calculations have 
been found in excellent agreement with experimental data, both for stable nuclei and for 
nuclei far away from the line of /3-stability. The NL3 interaction will also be used in the 
present analysis of ground state neutron densities. For the pairing field we employ the 
pairing part of the Gogny interaction 

1=1,2 

with the set DIS [|1^ for the parameters /Xj, Wi, Bi, Hi and Mj (i = 1,2). This force has 
been carefully adjusted to the pairing properties of finite nuclei all over the periodic table. 
In particular, the finite range of the Gogny force automatically guarantees a proper cut-off 
in momentum space. On the phenomenological level, the fact that we are using a non- 
relativistic interaction in the pairing channel of a relativistic Hartree-Bogoliubov model, 
has no influence on the calculated ground state properties. A detailed discussion of this 
approximation can be found, for instance, in Ref. [Q. 

The ground state of a nucleus results from a self-consistent solution of the RHB single- 
quasiparticle equations (|l|). The iteration procedure is performed in the quasi-particle basis. 
A simple blocking prescription is used in the calculation of odd-proton and/or odd-neutron 
systems. The resulting eigenspectrum is transformed into the canonical basis of single- 
particle states, in which the RHB ground state takes the BCS form. The transformation 
determines the energies and occupation probabilities of the canonical states. 

In Ref. we have performed a detailed analysis of ground state properties of Ni (28 < 
N < 50) and Sn (50 < N < 82) nuclei in the framework of the RHB model. In a comparison 
with available experimental data, we have shown that the NL3 + Gogny DIS effective 
interaction provides an excellent description of binding energies, neutron separation energies, 
and proton and neutron rms radii, both for even and odd-A isotopes. The RHB model 
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predicts a reduction of the spin-orbit potential with the increase of the number of neutrons. 
The resulting energy splittings between spin-orbit partners have been discussed, as well as 
pairing properties calculated with the finite range effective interaction in the pp channel. In 
Figs. |1| and ^ we plot the self-consistent ground state neutron densities for the even-A Ni 
(30 < N < 48) and Sn (56 < N < 74) isotopes. The density profiles display pronounced shell 
effects and a gradual increase of neutron radii. In the inserts we include the corresponding 
differences between neutron and proton rms radii. For Ni, the value of r„ — increases 
from ^ 0.03 fm for ^^Ni to 0.39 fm for ™Ni. The neutron skin is less pronounced for the Sn 
isotopes: r„ — rp = 0.27 fm for ^^^Sn. In Fig. ^ the difference of the neutron and proton rms 
radii for Sn isotopes are compared with recent experimental data [|1^]. The experimental 
values result from measured cross sections of the isovector spin-dipole resonances in Sn 
nuclei. The agreement between theoretical and experimental values is very good. The 
calculated differences are slightly larger than the measured values, but they are still within 
the experimental error bars. This result, together with the analysis of Ref. |0, indicates that 
the ground state neutron densities of the Ni and Sn isotope chains are correctly described 
by the RHB model with the NL3 + Gogny DIS effective interaction. 

The first experimental evidence of the formation of the neutron skin along a chain of 



stable and unstable isotopes was reported in Ref. [|13] for the Na nuclei. By combining 
rms charge radii, determined from isotope-shift data, with interaction cross sections of Na 
isotopes on a carbon target, it has been shown that the neutron skin monotonically increases 
to 0.4 fm in neutron-rich /3-unstable Na nuclei. In Ref. we have applied the RHB model 
in the description of properties of light nuclei with large neutron excess. The ground state 
properties of a number of neutron-rich nuclei have been analyzed: the location of the neutron 
drip-line, the reduction of the spin-orbit interaction, rms radii, changes in surface properties, 
and the formation of the neutron skin and the neutron halo. In particular, we have also 
calculated the chain of Na isotopes. The neutron density profiles are plotted in Fig. ^ with 
the differences of the neutron and proton rms radii in the insert. The calculations have been 
performed assuming spherical symmetry, i.e. deformations of proton and neutron densities 



were not taken into account. The blocking procedure has been used both for protons and 
neutrons. Strong shell effects are observed in the interior. In the central region the neutron 
density increases from ^ 0.065 fm"'^ for ^^Na, to more than 0.09 fm^^ for ^''Na. This effect, 
of course, corresponds to the filling of the s| neutron orbital. At r ^ 2 fm the density 
decreases from 0.1 fm~^ to 0.09 fm~^. The calculated neutron radii are compared with the 
experimental values fl^ in Fig. ^. The model reproduces the trend of the experimental data 



and, with a possible exceptions of = 11, the theoretical values are in excellent agreement 
with the measured radii. The smooth increase of r„ — Vp after = 11 (Fig. ^ is interpreted 
as the formation of the neutron skin. For A^ = 11 the number of protons and neutrons is the 
same, and therefore both the odd proton and odd neutron occupy the s — d orbitals with the 
same probabilities. An explicit proton-neutron short-range interaction, not included in our 
model, is probably responsible for a possible reduction of the neutron radius in this nucleus. 

In some loosely bound systems at the drip-lines, the neutron density distribution displays 
an extremely long tail: the neutron halo. The resulting large interaction cross sections have 



provided the first experimental evidence for halo nuclei The neutron halo phenomenon 



has been studied with a variety of theoretical models p9|-pT|. For very light nuclei in 
particular, models based on the separation into core plus valence space nucleons (three-body 
Borromean systems) have been employed. In Ref. the RHB model has been applied 
in the study of the formation of the neutron halo in the mass region above the s-d shell. 
Model calculations with the NL3 + Gogny DIS effective interaction predict the occurrence of 
neutron halo in heavier Ne isotopes. The study has shown that, in a mean-field description, 
the neutron halo and the stability against nucleon emission can only be explained with 
the inclusion of pairing correlations. Both the properties of single-particle states near the 
neutron Fermi level, and the pairing interaction, are important for the formation of the 
neutron halo. In Fig. ^ we plot the proton and neutron density distributions for ^°Ne, ^^Ne 
and ^^Ne. The proton density profiles do not change with the number of neutrons. The 
neutron density distributions display an abrupt change between ^^Ne and ^^Ne. A long tail 
emerges, revealing the formation of a multi-particle halo. The formation of the neutron halo 
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is related to the quasi-degeneracy of the triplet of states If 7/2, 2p3/2 and 2pi/2- For N < 22 
the triplet of states is in the continuum, it approaches zero energy at = 22, and a gap is 
formed between these states and all other states in the continuum. The pairing interaction 
promotes neutrons from the lf7/2 orbital to the 2p levels. Since these levels are so close in 
energy, the total binding energy does not change significantly. Due to their small centrifugal 
barrier, the 2p3/2 and 2pi/2 neutron orbitals form the halo. 

III. PARITY VIOLATING ELASTIC ELECTRON SCATTERING 

In this section we will illustrate how the neutron density distributions shown in Figs. 
|l] - P can be measured with polarized elastic electron scattering. We will calculate the 
parity-violating asymmetry parameter, defined as the difference between cross sections for 
the scattering of right- and left-handed longitudinally polarized electrons. This difference 
arises from the interference of one-photon and exchange. As it has been shown in Ref. 
0, the asymmetry in the parity violating elastic polarized electron scattering represents an 
almost direct measurement of the Fourier transform of the neutron density. 

The calculation procedure closely follows the derivation and definitions of Ref. 0. We 
consider the elastic electron scattering on a spin-zero nucleus, i.e. on the potential 

Vir)=V{r)+-f,A{r), (6) 

where V(r) is the Coulomb potential, and A(r) results from the weak neutral current ampli- 
tude 

The weak charge density is defined 

pwir) = J £r'GE{\r - r'|)[-p„(r') + (1 - Asm^ew)Pp{r% (8) 

where p„ and pp are point neutron and proton densities and the electric form factor of the 
proton is Ge{t) ^ ■^e"^'' with A = 4.27 fm~^. sin^GvK = 0.23 for the Weinberg angle. 
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In the hmit of vanishing electron mass, the electron spinor \& defines the helicity states 

^± = 1(1 ±75)^, (9) 

which satisfy the Dirac equation 

[a-p + y±(r)]^± = E^±, (10) 

with 

V±{r) = V{r)±A{r). (11) 

The parity- violating asymmetry Ai, or helicity asymmetry, is defined 

da+/dQ — dcr^/dQ , . 

' ^ da^/dVt + da./dVt' ^ ' 

where +(— ) refers to the elastic scattering on the potential V±{r). 

The calculation starts with the self-consistent relativistic Hartree-Bogoliubov ground 
state proton and neutron densities. The charge and weak densities are calculated by folding 
the point proton and neutron densities (see Eq. (H)). The resulting Coulomb potential 
V{r) and weak potential A{r) are used to construct V±{r). The cross sections for 
elastic electron scattering are obtained by summing up the phase shifts which result from 
the numerical solution of the partial wave Dirac equation. The calculation includes the 
Coulomb distortion effects. The cross sections for positive and negative helicity electron 
states are calculated, and the resulting asymmetry parameter Ai is plotted as a function of 
the scattering angle 9, or the momentum transfer q. 

In order to check the correctness and accuracy of the computer code which calculates 
the elastic electron scattering cross sections, we have performed the same tests as those 
reported in Ref. Q]: experimental data are reproduced for the elastic cross sections from 
^°^Pb at 502 MeV plane wave approximation results are reproduced; it is verified that 
the asymmetry parameter Ai is linear in the potential A{r) (|^. We have also reproduced 
the parity- violating asymmetries calculated in Ref. 0: elastic scattering at 850 MeV on 
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i6q^ ^^Ca, and ^ospi^ (relativistic mean-field densities), as well as on the three-parameter 
Fermi densities. 

In Figs. and ^ we plot the parity-violating asymmetry parameters Ai for the ^'^~'^^Ni 
isotopes, for the elastic electron scattering at 500 MeV and 850 MeV, respectively. The 
ground state neutron densities for these nuclei are shown in Fig. |1|. For electron energies 
below 500 MeV the asymmetry parameters are small (< 10^^), and the differences between 
neighboring isotopes are < 10~^. At 850 MeV (the energy for which most of the calculations 
of Ref. [§] have been performed), the values of Ai are of the order of 10~^. The differences 
between neighboring isotopes (~ 10~^) are especially pronounced for ^^~^^Ni, at ^ 20° 
and 6 ~ 35°. These differences reflect the strong shell effects calculated for the ground state 
neutron densities of the lighter Ni isotopes (see Fig. |l|). The heavier Ni isotopes display 
more uniform neutron densities in the interior, and the resulting asymmetry parameters Ai 
are not very different. Similar results are also obtained at 1000 MeV electron energy. Of 
course, above 1 GeV the approximation of elastic scattering on continuous charge and weak 
densities is not valid any more, and the structure of individual nucleons becomes important. 
In the remainder of this section we show the results for Ai at 850 MeV electron energy. For 
all chains of isotopes we have also calculated the asymmetry parameters at 250 MeV, 500 
MeV and 1000 MeV. The energy dependence, however, is similar to that observed for the 
Ni isotopes: below 850 MeV the differences in the calculated Ai for the neighboring isotopes 
are too small. 

The figures of merit 

for 850 MeV electron scattering on the Ni isotopes are shown in Fig. ^ The figure of 
merit is, of course, strongly peaked at forward angles (see also Fig. 12 of Ref. [Q]). In 
order to emphasize the differences between Ni isotopes, the F's are plotted in the interval 
10° < 9 < 30°. The figure of merit defines the optimal kinematics for an experiment in which 
the neutron density distribution could be determined from the measured parity-violating 

11 



asymmetries. 

The asymmetry parameter Ai provides a direct measurement of the Fourier transform 



of the neutron density M. In Figs. 10 and 11 we plot the asymmetries for the Ni isotopes 



(Fig. ^ as functions of the momentum transfer q = 2Esin6/2, and compare them with the 
squares of the Fourier transforms of the neutron densities 



4:71 r 

F{q) = — dr r^jo{qr)pnir). (14) 



The differences between the asymmetries can be directly related to the form factors. Note 
that the positions of the minima of Ai correspond almost exactly to the minima of the 
form factors. Of course, the agreement would have been perfect if we had plotted the 
Fourier transforms of the weak density (0), but the differences are indeed very small. More 
important is the observation that the measurement of the parity-violating asymmetry at 
high momentum transfer might provide information about the details of the density profile 
of the neutron distribution. 

In Fig. ^ we have shown that the differences in the neutron and proton rms radii, 
calculated with the RHB NL3 + Gogny DIS model, are in very good agreement with recent 
experimental data on Sn isotopes. In Fig. ^ we plot the asymmetry parameters Ai for the 
850 MeV elastic electron scattering on the even ^°^~^^^Sn isotopes. The angular dependence 
is similar to that observed for the Ni nuclei; significant differences between neighboring 
isotopes are only found at 9 > 20°. The asymmetry parameters are compared to the Fourier 



transforms of the Sn neutron densities in Figs. and [14 . 

The Na isotopes (Figs. ^ and display the formation of the neutron skin, as well as 
strong shell effects in the central region of neutron densities. These shell effects are clearly 
seen in the plots of the asymmetry parameters Ai in Fig. |T5[ Especially pronounced is the 
transition between ^^Na and ^''Na, which corresponds to the filling of the s| neutron orbital. 
The Fourier transforms of neutron densities are compared to the asymmetry parameters in 
Figs. and [T^. The differences between the minima of the asymmetry parameters for 
neighboring isotopes are the same as the differences between the minima of the Fourier 
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transforms of the densities. In principle, it should be possible to deduce the neutron density 
distribution from the measured asymmetries. 

An interesting theoretical question is whether the asymmetry parameters are sensitive 
to the formation of the neutron halo, i.e. whether parity violating electron scattering could 
be used to detect the formation of the halo. For the even Ne nuclei, the neutron halo 
phenomenon is illustrated in Fig. ^ The tail in the neutron density develops in ^^Ne. 
The Fourier transforms of neutron densities and the asymmetry parameters Ai are shown 
in Fig. At small momentum transfer the differences in the asymmetry parameters are 
very small. Only at g ^ 2.5 fm~^ the differences are of the order of 10^^. This probably 
means that, even if it became possible to measure polarized elastic electron scattering on 
extremely neutron-rich nuclei, the parity-violating asymmetries would not be sensitive to 
the formation of the neutron halo. 

IV. CONCLUSIONS 

The relativistic mean-field theory has been used to study the parity violating elastic 
electron scattering on neutron-rich nuclei. The parity-violating asymmetry parameter, de- 
fined as the difference between cross sections for the scattering of right- and left-handed 
longitudinally polarized electrons, provides direct information about the neutron density 
distribution. 

The ground state neutron densities of neutron-rich Ne, Na, Ni and Sn have been cal- 
culated with the relativistic Hartree-Bogoliubov model. The NL3 effective interaction has 
been used for the mean-field Lagrangian, and pairing correlations have been described by the 
pairing part of the finite range Gogny interaction DIS. The NL3 -|- Gogny DIS interaction 
produces results in excellent agreement with experimental data, not only for spherical and 
deformed /5-stable nuclei, but also for nuclear systems with large isospin values on both sides 
of the valley of /3-stability. In the present work, in particular, the calculated neutron rms 
radii are shown to reproduce recent experimental data on Na and Sn isotopes. 
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Starting from the relativistic Hartree-Bogoliubov solutions for the self-consistent ground 
states, the charge and weak densities are calculated by folding the point proton and neutron 
densities. These densities define the Coulomb and weak potentials in the Dirac equation 
for the massless electron. The partial wave Dirac equation is solved with the inclusion of 
Coulomb distortion effects, and the cross sections for positive and negative helicity electron 
states are calculated. The parity-violating asymmetry parameters are plotted as functions 
of the scattering angle 9, or the momentum transfer g, and they are compared with the 
Fourier transforms of the neutron density distributions. 

We have compared the parity- violating asymmetry parameters for chains of neutron-rich 
isotopes. For low electron energies (< 500 MeV), the differences between neighboring iso- 
topes are very small. At 850 MeV, significant differences in the parity- violating asymmetries 
are found at ^ > 20°. They can be related to the differences in the neutron density distri- 
butions. In particular, if plotted as function of the momentum transfer the asymmetry 
parameter can be related to the Fourier transform of the neutron density. It has been shown 
that the parity violating elastic electron scattering is sensitive to the formation of the neu- 
tron skin in Na, Ni and Sn isotopes, and also to the shell effects of the neutron density 
distributions. On the other hand, from the example of neutron-rich Ne nuclei, it appears 
that the asymmetry parameters would not be sensitive to the formation of the neutron halo. 
We conclude that, if it became possible to measure parity violating elastic electron scatter- 
ing on neutron-rich nuclei, the asymmetry parameters would provide detailed information 
on neutron density distributions, neutron radii, and differences between neutron and charge 
radii. This knowledge is, of course, essential for constraining the isovector channel of effec- 
tive interactions in nuclei, and therefore for our understanding of the structure of nuclear 
systems far from the /3-stability line. 
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Figure Captions 

FIG. 1. Self-consistent RHB single-neutron density distributions for even-N Ni (30 < N < 48) 
nuclei, calculated with the NL3 + Gogny DIS effective interaction. The differences between neutron 
and proton rms radii are shown in the insert. 

FIG. 2. Same as in Fig. |l[ but for Sn isotopes (56 < iV < 74). 

FIG. 3. Calculated differences between neutron and proton rms radii in Sn isotopes (squares), 
compared with experimental data from Ref. [13]. 

FIG. 4. Self-consistent RHB single-neutron density distributions for Na (12 < N < 21) isotopes. 

FIG. 5. Comparison of the calculated neutron rms radii for Na isotopes with experimental data 
from Ref. [14] . 

FIG. 6. Self-consistent RHB proton and neutron densities for 30,32,34]\j-g^ 

FIG. 7. Parity-violating asymmetry parameters Ai for elastic scattering from ^^"'^^Ni at 500 
MeV, as functions of the scattering angle 9. 

FIG. 8. Parity-violating asymmetry parameters Ai for elastic scattering from ^^^'^^Ni at 850 
MeV, as functions of the scattering angle 9. 

FIG. 9. Figures of merit ( p!3| ) as functions of the scattering angle 6, for elastic scattering from 
ss-'^^Ni at 850 MeV. 

FIG. 10. Parity-violating asymmetry parameters Ai (upper panel) and squares of normalized 
Fourier transforms of neutron densities (lower panel), as functions of the momentum transfer q, for 
elastic scattering from ^'^"Sep^j g^Q yi^y^ 
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FIG. 11. Same as in Fig. lid, but for elastic scattering from ^^Ni at 850 MeV. 



FIG. 12. Parity-violating asymmetry parameters Ai for elastic scattering from ^^^Sn at 850 
MeV, as functions of the scattering angle 9. 

FIG. 13. Parity-violating asymmetry parameters Ai (upper panel) and squares of normalized 
Fourier transforms of neutron densities (lower panel), as functions of the momentum transfer q, for 
elastic scattering from ^'^^"^-'^'^Sn at 850 MeV. 



FIG. 14. Same as in Fig. but for elastic scattering from ""^^^Sn at 850 MeV. 



FIG. 15. Parity- violating asymmetry parameters Ai for elastic scattering from ■^'^ ^■^Na at 850 
MeV, as functions of the scattering angle 9. 

FIG. 16. Parity-violating asymmetry parameters Ai (upper panel) and squares of normalized 
Fourier transforms of neutron densities (lower panel), as functions of the momentum transfer q, for 
elastic scattering from ^'^~^^Na at 850 MeV. 



FIG. 17. Same as in Fig. 16, but for elastic scattering from Na at 850 MeV. 



FIG. 18. Parity-violating asymmetry parameters Ai (upper panel) and squares of normalized 
Fourier transforms of neutron densities (lower panel), as functions of the momentum transfer q, for 
elastic scattering from 30,32,34^^^ g^Q j^gy^ 
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